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SUMMARY 

An analytical study has been made to obtain an idea of the lengths- 
needed for gaseous mixing in rocket combustion. Simplified models were 
assumed in which one gaseous propellant diffuses from several point 
sources into a uniform stream of the second gaseous propellant. No ef- 
fects of the combustion process were considered in these models. Degrees 
of mixing were calculated at several positions downstream by assuming 
homogeneous turbulence in the main stream. 

Detailed calculations were made for two different models. One assumes 
all the sources to be distributed in a single injection face. The other 
assumes sources distributed up to 7' inches downstream to simulate vaporiza- 
tion downstream of the initial injection point. 

For the first model good mixing is achieved for all the assumed tur- 
bulence conditions at a low ratio of downstream distance to source 
spacing; the estimated mixing lengths are less than lengths required for 
evaporation of typical rocket-engine sprays. Proper mixing takes signif- 
icantly longer for the downstream sources model. 

Calculations show that the contribution of molecular diffusion to the 
over-all mixing is negligible in these simple models, and mixing is there- 
fore controlled entirely by the turbulence parameters . This fact makes 
the obtained results independent of any particular propellant combination. 

These approximate calculations indicate that gaseous mixing generally 
occurs faster than spray evaporation. Its effect on the combustion rate 
certainly cannot be neglected, however, especially when sources are dis- 
tributed downstream of the first injection point. 
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OTTRODU CT ION 

The many physical processes which occur simultaneously in rocket com- 
bustion make the entire process extremely complex (ref. l) . The study 
of each process individually is one approach to identifying those that 
limit the rate of the combustion. The combustion chamber length required 
for evaporation of liquid-propellant sprays has already been studied 
(refs. 2, 3 , and 4). The present work estimates the combustion-chamber 
lengths required to mix propellants by turbulent and molecular gaseous - 
diffusion processes. 

A model and theory have been developed for calculation of gas con- 
centration profiles for systems in which both propellant vaporization and 
gaseous diffusion occur (ref. 5, p. 76). To satisfy the purposes of this 
preliminary survey, the complete treatment of reference 5 has not been 
used. Instead calculations for turbulent and molecular diffusion of gases 
have been carried out for simplified analytical models. The results are 
then compared with those from the vaporization calculations. 
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SYMBOLS 

concentration of diffusing gas, Ib/cu ft 
turbulent diffusion coefficient, sq ft/sec 
molecular diffusion coefficient, sq ft/sec 
scale of turbulence, ft 
ambient pressure, Ib/sq in. abs 

radial distance in a y-z plane from any point in the stream to the 
axial line through a source 

source spacing in the y-z plane, in. 

turbulence intensity, V^2/U 
constant mean stream velocity, ft/sec 

root-mean-square turbulent -velocity fluctuation in any direction 
source strength, lb/sec 


x 


distance downstream of injection plane 
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a T(x/s) 

a standard geometric deviation 

o 

\|r ripple factor 

0 ) spreading factor, sq ft 

Subscripts : 
max maximum 

min minimum 


THEORY AND MODELS FOR THE CALCULATIONS 
General Theoretical Equations 

The differential equation for combined molecular and turbulent mass 
transfer in three dimensions is (ref. 5, p. 63): 

Sc _ dfc dfc c^C , ^ 

Sx^ 

where C is the concentration of the diffusing substance at a given posi- 
tion and time, and a), the time-dependent spreading coefficient, is de- 
fined by 

CJO = f (D + d) dt* (2) 

* / 0 


In the models to be described later all diffusion is assumed to come from 
a number of continuous point sources. Each point source is assumed to 
contribute independently to the total concentration of the diffusing gas 
at any given point. A good approximate solution of equation (l) for a 
continuous point source is given in reference 5 as 


C = _L_ e -r 2 /4co 
W7U 4nm 


( 3 ) 


This equation pertains to the diffusion model shown in the following 
sketch: 
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Material diffuses into an infinite stream of uniform density flowing in 
the x-direction with linear velocity U. The equation involves three 
assumptions: (l) the source material occupies negligible volume in the 

main stream, (2) diffusion in the x-direction is negligible, and (3) r is 
very much smaller than x. In connection with the last assumption, how- 
ever, inspection of the data of reference 6 (p. 16) indicates that the 
error involved in using this equation at r/x values up to 0.4 is negli- 
gible for the purposes of this report. The spreading factor cjo is evalu- 
ated at time t = x/U. A simple analytical expression for cd obtained 
under the assumption of homogeneous turbulence is (ref. 5, p. 63) 

CD = LxT + L 2 (e" xT / L - l) + (4) 

where L is the scale of the turbulence and T is the turbulence inten- 
sity, defined as the ratio of the root-mean-square velocity fluctuation 
to the mean stream velocity. 

It should be noted that the concentration C in equation (3) is a 
time average value and tells nothing about the fluctuation of concentra- 
tion with time at any point. It should also be noted that the effect on 
mixing of the combustion process is not being considered in these 
equations . 
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For diffusion from a number of point sources, mixing at any point 
can most readily be shown in terms of a Tt ripple factor.” This factor is 
defined as 


♦ = c — 

°max 

The terms C m i n and C max are the lowest and highest values of the 
diffusing-material concentration at the particular downstream distance. 

The ripple factor for perfect mixing is unity. 

It can be pointed out here that perfect mixing is unnecessary for com- 
plete combustion except at stoichiometric propellant ratios. Mixing in 
a combustion process has to proceed only to a point where one propellant 
has been completely consumed. However, a high degree of mixing is impor- 
tant in achieving a uniform temperature profile in the hot gases, and 
thus, high performance. 


( 5 ) 


Models Used 

Two models were used in these calculations. They are illustrated in 
figures 1 and 2 in which some typical concentration profiles are shown 
qualitatively. These models are as follows: 

Model 1: An infinite rectangular grid of equal strength, evenly 

spaced point sources in the x = 0 plane. 

Model 2; A lattice made up of eight grids of sources as described 
for model 1. The grids are spaced at 1-inch intervals downstream along 
the axis of flow with the first one at the x = 0 plane . 

The time -average diffusing substance concentration at any point is 
assumed to be the sum of contributions from all sources acting independ- 
ently at that point, so that equation (3) can be used. The results apply 
only for the situation where no burning takes place and where the 
diffusing-material volume is small compared with the mainstream volume. 
All calculations performed with these models neglect any boundary effects. 
(These wall effects can be calculated using the theory of ref. 7.) 


METHOD OF CALCULATION 

Assignment of Turbulence Parameters 

The calculation of downstream concentration profiles using the pre- 
ceeding equations requires numbers for the turbulence scale and intensity. 
Here a temporary impasse appears, inasmuch as these values are not known 
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for the rocket-combustion environment. Preliminary calculations showed 
that turbulent mixing is not affected by the scale (the over-all size of 
the eddies) except at very small values. Therefore, a scale value of 0.1 
foot, based on pipe-turbulence data, was arbitrarily used throughout these 
calculations. A series of turbulence-intensity values was used. The 
intensities used, and also the downstream distances and source spacings 
(in the x = 0 plane), are as follows: 

T = 0.05, 0.07, 0.09, 0.11 
x = 1 to 12 in. at intervals of l/2 or 1 in. 
s = l/2 and also 1 in. for model 1 only 


Relative Contributions of Molecular and Turbulent Diffusion 

At this point it is possible to evaluate the relative contributions 
of turbulent and molecular diffusion to the entire mixing process by cal- 
culating the two separate contributions to the spreading factor <u. This 
task requires using a particular propellant combination and assuming 
specific flow conditions. The system oxygen diffusing into hydrogen has 
been used, because this combination has the largest molecular diffusion 
coefficient at any given temperature and pressure. Furthermore, to exag- 
gerate the contribution of the molecular diffusion term to the spreading 
factor 00 , low values of turbulence level were assumed as follows: 


Scale of turbulence, L, ft 0.10 

Turbulence intensity (normal pipe flow value), T 0.03 

Ambient pressure, p, lb/sq in. abs 600 

Constant mean stream velocity, u, ft/sec 100 

Average molecular-diffusion temperature, °R 400 


The molecular-diffusion coefficient for these conditions is 1.207 x 10“5 
square feet per second as calculated by the equation of reference 8. The 
result of this calculation shows that the maximum contribution of the 
molecular diffusion term to the whole spreading factor is only 0.3 percent 
of the turbulent -mixing contribution under the assumed conditions. It can 
be concluded that molecular diffusion is always negligible for the models 
employed, regardless of the propellant combination chosen. 
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Calculations Involving the Ripple Factor 


Analytical expressions for the ripple factor can be obtained as a 
result of neglecting molecular diffusion. The general equation for the 
ripple factor for n continuous point sources is 


^min 

^max 



( 6 ) 


where 


r = r value for C min 


a = number of sources having a given value of r 


r' = r value for 

b = number of sources having a given value of r T 

The equation can be simplified in its application to model 1 by neglecting 
all but the most important source contributions . Details are given in the 
appendix for the derivation of an approximate equation for model 1. The 
final result and the assumptions used to obtain it are as follows: 

Assumption 1: Only sources with r = s^Jz/z, s^fTti/z, and s^/18/2 
contribute to C m -^ n . 

Assumption 2: Only sources with r f =0, s, and contribute to 

Cmax • 

The resulting equation is 


4e -(l/4a 2 ) + 8 


5 

+ 4 1 

[ e -(l/ta 2 )] 9 

1 + 4 

' e -(l/4a 2 )J 

r-« 

| e -(l/4a 2 )] 

I 4 


where a = T(x/s) 
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RESULTS AND DISCUSSION 
Model 1 

Figures 3(a) and (b) show results of the detailed ripple-factor cal- 
culations for this grid of sources model. At the l/2-inch source spacing 
(fig. 3(a)) the ripple factor becomes equal to 0.95 at an x value of 
5.8 inches for the lowest turbulence intensity used (T = 0.05). This 
fact indicates that mixing due to turbulence is fairly rapid for a source 
spacing typical of injector port spacings in a rocket engine. 

Mixing occurs much more slowly for the 1-inch spacing (fig. 3(b)). 
This difference in the degree of mixing is resolved in figure 4 where the 
ripple factor is plotted as a function of the ratio of downstream distance 
to source spacing. For a fixed turbulence intensity these plots show that 
the degree of mixing is a function of the ratio x/s, within an accuracy 
of 10 percent . 

The dependence of ripple factor on turbulence intensity is shown in 
figure 5. Obviously mixing increases with increasing intensity. 


Comparison with Approximate Equation 

The approximate equation for ripple factor (eq. (7)) predicts that 
mixing is independent of turbulence scale for this model and depends 
only on the turbulence intensity and. the ratio of downstream distance to 
source spacing (x/s). This equation is shown in the semilogarithmic plot 
of figure 6 with the curve for the detailed calculations. In these plots 
high values of a = T(x/s) give good mixing. The comparison shows that 
equation (7) does a good qualitative job of describing the mixing over a 
wide range of turbulence conditions. In general, the quantitative agree- 
ment is not good, but improves as the ripple factor becomes greater than 
0.7. The equation is therefore very useful in predicting quickly the 
point where complete mixing is approached. 


Model 2 

Ripple factors for the three dimensional lattice of point sources are 
shown in figure 7. This model approximates the diffusion of a liquid 
which is vaporizing uniformly as it travels downstream. (As mentioned in 
the INTRODUCTION an exact treatment of this type of calculation is given 
in ref. 5.) At a turbulence intensity of 0.05 with l/2-inch spacing, a 
downstream distance of 11.3 inches is now required to obtain 95 percent 
of perfect mixing, that is, 4 inches beyond the last downstream sources. 
This compares with the 5.8 inches needed for model 1. As long as sources 
are present downstream, the ripple factor increases rather slowly. 
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However, downstream of the last sources the ripple factor reaches unity 
in a shorter distance for the lattice model than it does for the grid of 
sources. This is shown by the curves in figures 4 and 7. 


Comparison of Mixing and Vaporization Length Calculations 

The approximate mixing lengths just calculated can be compared with 
the results of vaporization-rate calculations reported in references 3 
and 9. A comparison of the lengths required for 95-percent mixing with 
those required for 95-percent vaporization of the mass of typical rocket- 
engine sprays will give an idea of the relative importance of the two 
processes when they occur simultaneously. The table below summarizes cal- 
culations from references 3 and 9 for several heptane sprays and one oxygen 
spray evaporating into an assumed rocket-engine combustion environment. 


Mass 

median 

drop 

radius, 

in. 

Standard 

geometric 

deviation, 

a g 

Length for 95- 
percent mass 
vaporized, 
in. 

Heptane 

Oxygen 

0.003 

1.25 

5.9 




1.54 

7.8 




2.3 

14 

5.8 



3.6 

25 



The larger the value of cjg, the greater is the radius of the largest drop 
in the spray for a given mass median drop radius. The largest 95-percent 
mixing lengths found in the present calculations (at the l/2-inch source 
spacing) are 5.8 inches for model 1 and 11.3 inches for model 2. This 
information indicates that for cases with Cg >1.25, gaseous mixing from 

a single grid of point sources is faster than vaporization of a heptane 
spray. Reference 3 shows that several actual injector sprays have a cjg 
value of about 2.5. The conclusion is that gaseous mixing for the grid 
of sources model occurs in a shorter length than is needed to vaporize 
typical rocket-engine heptane sprays. Additional calculations in refer- 
ence 9 have shown this statement to be true for fuels other than heptane. 
The vaporization length for the oxygen spray is the same as the maximum 
mixing length found for the grid of sources model. So if a fuel more 
volatile than oxygen is used with it, the mixing rate still does not be- 
come the limiting factor on the entire combustion rate, even at the lowest 
turbulence intensity used. 

The mixing length calculated for a lattice of diffusion sources ex- 
tending 7 inches downstream is not significantly less than the vaporization 
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length for a heptane spray with a a g value of 2.3. No quantitative 

significance is given to the actual values. The comparison does show that, 
although the mixing process is quite rapid downstream of the last diffusion 
sources, it is by no means fast enough to be completely neglected as a 
combustion-rate controlling factor. 

The preceeding comparison of the rates of mixing and vaporization has 
been made on the basis of time-average concentration calculations. The 
time variation of concentration at a given point can be quite significant. 
(The root -mean-square value of the fluctuation can be as high as 80 per- 
cent of the mean value.) What effect these time fluctuations will have 
upon the mixing process is not known from the present work. (See ref. 

10 for a method of treating these fluctuations.) 

The several principle assumptions used in this analysis and their 
implications with respect to a rocket-combustion process require some 
further discussion: 

(1) The volume occupied by the diffusing gas is assumed negligible 
compared with that of the main stream. In a practical case this assumption 
is, of course, never realized. For example, in the combustion of hydro- 
gen and oxygen at an oxidant-fuel mass ratio of 4.0 the ratio of oxygen 

to hydrogen volumes is 1:4. In a constant area duct, therefore, the in- 
jection of gaseous oxygen would accelerate the total flow by approximately 
25 percent . 

(2) The mixing is assumed to occur between nonreacting gases. The 
contribution of the combustion process to the mixing is not considered. 

The idea of "f lame -generated" turbulence was proposed several years ago, 
but no experimental evidence has yet been found to support the idea that 
combustion significantly increases the scale and intensity of turbulent 
flow (see refs. 11 and 12). The fact remains, however, that little is 
known about the effect of the combustion process on mixing. 

(3) The turbulence scale and intensity are assumed to remain constant 
throughout the mixing process. This assumption is a direct consequence 

of assumption (2) . Photographic studies of combustion indicate that flow 
downstream of the principal combustion zone is pseudo laminar (ref. 13). 

It seems likely that, in a real combustor, turbulence intensity diminishes 
rapidly beyond the main reaction zone. 


SUMMARY OF RESULTS 

For gaseous diffusion from groups of continuous point sources the 
following results are summarized: 
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1. Turbulent mixing is practically independent of the scale of turbu- 
lence except at very low values. 

2. Mixing increases with increasing turbulence intensity and with in- 
creasing ratio of downstream distance to source spacing. 

3. On the basis of these time-average concentration calculations, tur- 
bulent mixing is faster than the vaporization process for the simple grid 
of sources model at all the turbulence intensities used. This work has 
told nothing, however, about the time -dependent concentration fluctuations, 
which can be large enough to affect the mixing process. 

4. Practically all the mixing occurs by turbulent diffusion in the 
models used. Molecular diffusion contributes negligibly to the mixing 

so that the results of the present work are independent of any particular 
propellant combination. This is true for source types other than the 
continuous point source. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, June 10, 1958 
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APPENDIX - ANALYTICAL EXPRESSION FOR THE RIPPLE FACTOR IN MODEL 1 

Figure l(b) shows the geometric relations for this model, as viewed 
along the x-axis. Only points near the center of the grid are considered. 
Point A is a typical point on a source streamline for the calculation of 
Cmax* The most important sources contributing at this point are (l) one 
with r f = 0, (2) four with r T = s, and (3) four with r T = s^2. The 
typical point for C min is point B. In this case the contributing sources 

are (l) four with r = s^fz/2 (indicated by the squares), (2) eight sources 
with r = s^/lO/2 (indicated by the triangles), and (3) four with 
r = s^/l8/2 (indicated by the diamonds) . If only these three values of 
r and r T are used in the general equation (6), it becomes 

ie- 52 / 8 ^ + 8e- 5s2 / to + 4e- 9s2 / 8Q ( ... 

1 ♦ ♦ «.-*/»> 

Next, equation (4) for cd can be simplified by dropping the molecular 
diffusion term, expanding the exponential, and dropping all terms above 
second order. The result is 


CD 


x 2.ji2 

2 


Substitution of equation (A2) into equation (Al) gives 

+ 8 ( e -l/ta 2 ) 5 + 4 ( e -l/ta 2 ) 9 

1 + 4 ( e -^ 2 ) 2 + t ( e - V ^ 2 ) 3 4 


(A2) 


(7) 
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Figure 1. - Diffusion model for a grid of evenly spaced point sources, model 1. 
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Figure 2. - Diffusion model for a lattice of point sources, model 2. 
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(a) Source spacing; l/2 inch. 

Figure 3. - Ripple factor against downstream distance for the grid 
of sources, model 1. Turbulence scale, 0.10 foot. 
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(b) Source spacing, 1 inch. 

Figure 3. - Concluded. Ripple factor against downstream distance 
for the grid of sources, model 1. Turbulence scale, 0.10 foot. 
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Figure 4. - Ripple factor against distance to spacing ratio for 
model 1. 
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Figure 5. - Ripple factor against turbulence intensity for the grid 
of sources, model 1. 
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Figure 6. - Approximate ripple factor equation for model 1. 
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Figure 7. - Ripple factor against distance to spacing ratio 
for model 2 } source spacing (in y-z plane), l/2 inch. 
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